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ABSTRACT 

We construct Super-Eddington Slim Disks models around both stellar and super- 
massive black holes by allowing the formation of a porous layer with a reduced effective 
opacity. We show that at high accretion rates, the inner part of the disks become 
radiation pressure dominated. However, unlike the standard scenario in which the 
disks become thick, super-Eddington disks remain slim. In addition, they accelerate 
a significant wind with a "thick disk" geometry. We show that above about 1.5 times 
the standard critical mass accretion rate (needed to release the Eddington luminosity) , 
the net luminosity released is above Eddington. At above about 5 times the standard 
critical rate, the central BH accretes more than the Eddington accretion rate. Above 
about 20m crit , the disk remains slim but the wind becomes spherical, and the present 
model breaks down. 

Key words: Accretion disks 



1 INTRODUCTION 

The accretion of matter onto compact objects can often be 
described using the run-of-the-mill thin accretion disk model 
of Shakura & Syunyaev (1973, S&S). Because the accretion 
disk is optically thick, matter can radiate the potential en- 
ergy it dissipates and remain cold, thus forming a geomet- 
riclaly "thin" disk. The turbulent viscosity responsible for 
the dissipation is often described through the standard a- 
model, and it is also responsible for the transport of angular 
momentum outwards. 

The S&S thin disk model applies to a wide range of con- 
ditions found in nature, however, when one of the underlying 
assumptions break down, so does the model. 

At sufficiently low accretion rates, the disk becomes 
optically thin and it cannot radiate the energy dissipated. 
This energy is therefore advected with the flow, forming the 
so called Advection Dominated Accretion Flow (Ichimaru 
1977; Narayan & Yi 1994; Abramowicz et al. 1995). Because 
of the high temperatures and pressures, ADAFs inflate to 
become geometrically "thick" and sub-Keplerian. Another 
interesting aspect of disks in this accretion regime is the 
possibility of generating significant outflows (Blandford & 
Begelman 1999; Shaviv et al. 1999). 

The inability to radiate enough energy also arises for 
very high accretion rates, giving rise to advection dominated 
flows once the disk becomes radiation pressure dominated 
(Paczynsky & Wiita 1980; Jaroszynski et al. 1980). 

The energy release rate in a Keplerian disk, down to a 
radius r, is given by GMB.w_m/2r, where Mbh is the black 
hole mass, and m is the mass accretion rate. This energy 



should be compared to the Eddington luminosity of the BH, 
defined as: 



JEdd 



47TCGM B H 



(1) 



from which a critical accretion rate can be define to be 

-t/Edd 



W-crit = 



(2) 



with rjo — 1/16 being the standard efficiency for accretion 
around a Schwarzschild BH. If m is large enough, that is, 
m > m cr it then the energy released will approach the Ed- 
dington luminosity before reaching the inner radius of the 
disk. From that point inwards, the radiation pressure will be 
dominant. The disk will become geometrically thick (having 
a scale height H(r) ~ r), and the local flux will approach 
the Eddington flux defined as 



cGMbhz 
kR 3 ' 



(3) 



with R = \/r 2 + z 2 (where r and z are the cylindrical co- 
ordinates). However, because the Eddington flux cannot be 
surpassed, the dissipated energy cannot be entirely radiated 
away and part of it should therefore be advected inwards. 
The high pressure also implies that the disk rotates at sub- 
Keplerian velocities. 

Other interesting aspects of these disks is their possi- 
ble ability to accelerate very luminous jets (Abramowicz & 
Piran 1980), and the fact that although the local flux in 
the disk does not surpass the Eddington flux, the overall lu- 
minosity can surpass the Eddington luminosity, simply due 
to the disk geometry (Jaroszynski et al. 1980; Paczynsky & 
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Wiita 1980). Roughly, the Eddington luminosity can be sur- 
passed by a factor of ~ In r ou t/r in , where the radii denote 
the inner and outer extents of the radiation pressure disk. 

For nearly critical accretion, an intermediate type of 
solutions exists, that of "slim-disks". However, because it 
exists only near the last stable orbit (thus allowing for ad- 
vection of heat into the Roche Lobe overflow), its solution 
requires relativistic corrections (Abramowicz et al. 1988). 

The above models assume, however, that the Eddington 
flux cannot be surpassed locally. Nevertheless, it was shown 
that super-Eddinigton states do naturally arise in nature 
(Shaviv 2001b), allowing for high luminosities, while gener- 
ating optically thick winds. Our goal in the present work 
is to consider the recent advances in the understanding of 
how super-Eddington atmospheres arise, and what they look 
like, and to incorporate these ideas into models for very high 
accretion rate accretion disks. 

In §2, we begin by reviewing our present understanding 
of how super-Eddington states arise. In §3 we describe our 
model for super-Eddington accretion, and in §4 we describe 
the numerical solution. In §5 we describe the numerical re- 
sults, and end with a discussion in §6. 



2 BACKGROUND: SUPER-EDDINGTON 
STATES 

Wc begin by reviewing the relevant physics pertaining to the 
emergence of super-Eddington states. These include three 
particular elements. First, the rise of inhomogeneities due 
to radiative-hydrodynamic instabilities was shown to re- 
duce the effective opacity (Shaviv 1998). This allows for 
the existence of super-Eddington atmospheres. Second, once 
a super-Eddington state arises, strong continuum driven 
winds are accelerated (Shaviv 2000, 2001b; Owocki et al. 
2004). These optically thick winds are responsible for a sig- 
nificant mass loss and are also important when determining 
the appearance of these objects. Last, if the wind mass loss 
is too large, wind stagnation and a photon-tired state arises 
(Owocki & Gayley 1997; van Marie et al. 2009). In it, a 
layer is formed in which strong shocks mediate a high energy 
flux without an excessive mass flux. These three components 
are the necessary building blocks for the Super-Eddington 
(SED) accretion disk models, and we therefore review them 
below. 

2.1 The rise of super-Eddington states 

According to common wisdom, objects cannot shine beyond 
their classical Eddington limit, £Edd (or locally beyond the 
Eddington flux, J^dd, as is the case in accretion disks), since 
no hydrostatic solution exists. In other words, if objects do 
pass ^Edd, they are highly dynamic. They have no steady 
state, and a huge mass loss should occur since their atmo- 
spheres are then gravitationally unbound and they should 
therefore be expelled. Thus, astrophysical objects according 
to this picture, can pass £Edd but only for a short dura- 
tion corresponding to the time it takes them to dynamically 
stabilize once SED conditions arise. 

For example, this can be seen in detailed ID numeri- 
cal simulations of thermonuclear runaways in classical nova 
eruptions, which can achieve SED luminosities but only for 



several dynamical time scales (e.g., Starrfield 1989). How- 
ever, once they do stabilize, they are expected and indeed 
do reach in the simulations, a sub-Eddington state. Namely, 
we naively expect to find no steady state SED atmospheres. 
This, however, is not the case in nature, where nova erup- 
tions are clearly SED for durations which are orders of 
magnitude longer then their dynamical time scale (Shaviv 
2001b). This is exemplified with another clear SED object — 
the great eruption of the massive star 77-Carinae, which was 
a few times above Eddington for over 20 years (Shaviv 2000) 
The existence of a super-Eddington state can be natu- 
rally explained, once we consider the following: 

(i) Atmospheres become unstable as they approach the 
Eddington limit. In addition to instabilities that operate 
under various special conditions (e.g., Photon bubbles in 
strong magnetic fields, Arons 1992; Gammie 1998; Begel- 
man 2002, or s-mode instability under special opacity laws, 
Glatzel 1994; Papaloizou et al. 1997), two instabilities op- 
erate in Thomson scattering atmospheres (Shaviv 2001a). 
It implies that all atmospheres will become unstable already 
before reaching the Eddington limit. 

(ii) The effective opacity for calculating the radiative 
force on an inhomogeneous atmosphere is not necessarily 
the microscopic opacity. Instead, it is given by 



where ( ) v denotes volume averaging and F is the flux (Sha- 
viv 1998). The situation is very similar to the Rosseland vs. 
Force opacity means used in non-gray atmospheres, where 
the inhomogeneities are in frequency space as opposed to 
real space. For the special case of Thomson scattering, the 
effective opacity is always reduced. 

Thus, we find that as atmospheres approach their clas- 
sical Eddington limit, they will necessarily become inhomo- 
geneous. These inhomogeneities will necessarily reduce their 
effective opacity such that the effective Eddington limit will 
not be surpassed even though the luminosity can be super- 
classical-Eddington. This takes place in the external regions 
of luminous objects, where the radiation diffusion time scale 
is shorter than the dynamical time scale in the atmosphere. 
Further inside the atmosphere, convection is necessarily ex- 
cited such that the total energy flux may be SED, but the 
radiative part of it is necessarily sub-Eddington with the 
convective flux carrying the excess (Joss et al. 1973). 

2.2 Super-Eddington Winds 

The atmospheres of SED objects, and SED accretion disks 
in particular, effectively remain sub-Eddington while being 
classically super-Eddington, only as long as the inhomo- 
geneities comprising them are optically thick. This condition 
will break at some point where the density is low enough. 
At this height, the effective opacity returns to its micro- 
scopic value and hence the radiative force becomes super- 
Eddington again. From this point outwards we obtain con- 
tinuum driven winds. Because the winds are generally op- 
tically thick, the conditions in them affect the structure of 
the disk beneath. 

At the critical point, the radiative and gravitational 
forces balance each other. This point will coincide with a 
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sonic surface for a steady state wind (where the mass loss 
velocity equals the local speed of sound). This allows us to 
obtain the local mass loss rate per unit area, given by 



0wind = Pcrit«s(2atm) = COUSt., 



(5) 



where 2atm is the vertical height of the critical point, p C rit is 
the density at this point and v s is the local speed of sound. 

Based on the fact that instabilities develop structure 
with a typical size comparable to the density scale height in 
the atmosphere, it is possible to estimate the average density 
at the sonic point (Shaviv 2001b). Using this density, the 
mass loss can be estimated to be 

, F — JEdd 



ind 



= w- 



(6) 



where W is a dimensionless wind "function" . In principle, 
W can be calculated ab initio only after the nonlinear state 
of the inhomogeneities is understood. This however is still 
lacking as it requires elaborate 3D numerical simulations of 
the nonlinear steady state. 

Nevertheless, deriving W can be achieved in several phe- 
nomenological models which depend on geometrical param- 
eters such as the average size of the inhomogeneities in units 
of the scale height (£ = d/l p ), the average ratio between the 
surface area and volume of the blobs in units of the blob 
size (H), and the volume filling factor / of the dense blobs. 
For example, in the limit in which the blobs are optically 
thick, one can show that W w 3H/320;/C(l - f) 2 (Shaviv 
2001b), with v being the ratio between the effective speed 
of sound in the atmosphere to the adiabatic one. Thus, W 
depends only on geometrical factors. It does not depend ex- 
plicitly on the Eddington parameter Y = L/C^dd as long as 
the blobs have a single length scale. Once this assumption 
is alleviated, W can become a weak function of T (Owocki 
ct al. 2004). Comparison to observations yields typical val- 
ues of W ~ 1 - 10 (Shaviv 2001b). 

2.3 Photon-Tired Winds 

An interesting modification to the above continuum driven 
winds arises when the predicted mass loss is too high for 
the available luminosity to push it to r — > oo. This happens 
when v esc > yJv a c/W, and it gives rise to "photon-tired 
winds" (Owocki & Gayley 1997). A wind solution with a 
monotonically decreasing velocity is then not possible, be- 
cause the wind stagnates at a finite radius. 

The behaviour of photon tired winds was studied by 
van Marie et al. (2009). It was found that shocks form be- 
tween infalling material and the outflowing wind. This forms 
a layer of shocks in which there is a large kinetic flux, but 
without the associated mass flux. When photon tired winds 
arise, the mass loss from the top of the layer of shocks is re- 
duced to less than the photon tired limit, and the luminosity 
to less than the Eddington luminosity. 



3 THE MODEL 

Our goal is to construct Super-Eddington accretion disks, 
namely, disks which radiate fluxes that can locally exceed 
the Eddington flux. We look for slim disk solutions which 
are heuristically described in fig. 1. More specifically, we 
have to consider the following points. 



• Geometry: Although in principle, one could envision 
super-Eddington accretion solutions with different geome- 
tries, we look for disk like solutions. We shall assume that 
the vertical length scale is sufficiently smaller than the ra- 
dius, such that we can deconstruct the problem into vertical 
and radial components. As we shall see below, the hydro- 
static part of the disk satisfies this conditions except for 
the highest accretion rates. Note that because the vertical 
structure is not much smaller, the geometry is not that of a 
thin disk, but that of a slim disk. This also implies that we 
cannot assume Keplerian velocities. 

• Porosity: As mentioned above, high radiative fluxes give 
rise to porosity which reduces the effective opacity, thereby 
allowing the existence of super-Eddington fluxes in the hy- 
drostatic atmospheres. Thus, a necessary component of our 
model is an opacity law of the form k(T), which takes the 
porosity into account. 

• Convection: Joss et al. (1973) have shown that high 
radiative fluxes give rise to convection as the radiative fluxes 
approach the Eddington limit. The dense inner parts of the 
SED disk, near the equatorial plane are therefore expected 
to be convective. 

• Wind: Since a porous atmosphere can reduce the opac- 
ity only as long as the inhomogeneities comprising it arc 
optically thick, a necessary outcome of SED atmospheres is 
the acceleration of continuum driven winds where the av- 
erage density is low enough. Because this mass loss can be 
significant, it has two interesting ramifications. First, be- 
cause the wind is generally optically thick, the photosphere 
is going to reside in the wind. This will have various ob- 
servational consequences. Second, the mass accretion rate 
decreases as the radius decreases. 

3.1 Radial Structure 

As mentioned above, the small thickness of the disk enables 
the separation between the radial structure and the vertical 
one. The equations describing the radial structure are ob- 
tained from the radial conservation of mass, radial momen- 
tum, energy, and angular momentum. An additional equa- 
tion is the closure relation for the stress tensor. 
Radial mass conservation gives: 



dm 



= 47rr0 wi nd, 



dr 

where 

rh = Amrv r I pdz 



Jo 



(7) 



(8) 



is the mass accretion rate. Note that we assume a height 
independent velocity structure. 

Conservation of radial momentum gives 

v dvr + \dF_ _ 

dr p dr dr ' 

where * is a pseudo-Newtonian potential given by Paczyn- 
sky & Wiita (1980) as 

* = -f^, (10) 
R-rg 

and R = \fr 2 + z 2 . For the radial structure, we assume that 

2 = 0. 

The equation for angular momentum conservation is 
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Figure 1. The model structure. The different regions are: (A) A sub-Eddington thin disk (following the solution of S&S). (B) The 
accretion disk loses its flat geometry at the radius where the energy release corresponds to the Eddington luminosity. The disk inflates 
and becomes radiation-pressure dominated. Once super-Eddington states are allowed to arise, the standard thick disk picture is modified. 
(C) A porous layer forms at the less dense regions above the convection layer. (D) At a height were the porous structure become optically 
thin, a wind is accelerated. (E) Since it is optically thick, the photosphere is located in the wind itself. (F) In the inner parts of the disk, 
the escape velocity is large enough to give rise to a photon tired layer, which effectively moves the sonic point higher. 



where r rr p is the tangential stress. Following the standard 
"a" prescription, we write: 



= -aP. 



(12) 



The heat produced by the viscosity is partly radiated 
away (locally), and partly transferred by advection into 
smaller radii. The radial heat advection is given by the dif- 
ference between the amount of heat being produced by the 
viscosity, and the energy radiated away from the surface of 
the disk, that is, 



dr m 



(13) 



where s is the specific entropy, F is the flux radiated from 
the surface, while <£> is the dissipation function, given by 



$ — r 



du) 
dr 



r rtf) dz. 



(14) 



Note that the specific entropy increases with decreasing radii 
in the case of advection, hence ds/dr < 0. However, another 
possibility exists, in which the infalling matter releases part 
of its heat (i.e., ds/dr > 0). This heat, together with the heat 
generated by the viscosity, is radiated away. This process 
dominates the inner radii region of the accretion flow, as we 
will show in the results. 



3.1.1 Inner Sonic Point 

At large radii, the radial velocity is very small when com- 
pared to the speed of sound, while matter is freely infalling 
at the vicinity of the BH. This implies that at some radius 
r — r s , the radial velocity should become equal to the speed 
of sound, and then exceed it. Similarly to the case of spher- 
ical Bondi accretion, we can obtain the condition for this 
point. Because of the radial pressure gradient, this point 
will reside between the last stable orbit r; s = 3r 9 and the 
marginally stable orbit, r ms = 2r g . Using eqs. 8-9, and that 



dp dp ds dp dp 

dr ds dr dp dr s V3 J dr B dr 



(i->) 



ds 2 dp 

v. 



(15) 



where vl = dp /dp, and s is the specific entropy, we obtain 

vl - V 2 S dv r 1 r 2 , , 2 2 dp ds 1 

~' - [ V s + K v <j> ~ «kep)J 



dr 



ds dr p 



(16) 



Note that we have assumed the following approximation, 
that 



(17) 



1 dE _ 1 dp 
E dr p dr 

where E = J pdz. At the sonic radius v r = v s , such that the 
r.h.s. of eq. 16 must vanish, i.e., at the sonic radius we have 
that 

> ds 1 



« s + fa - v kcp ) ~ r 



ds dr p 



(18) 



We use this expression as the inner boundary condition for 
the disk. 



3.2 Vertical Structure 

The vertical structure of the SED accretion disk can be di- 
vided into two regions, a hydrostatic region which includes 
also the porous atmosphere, and the region of a continuum 
driven wind. Because the wind is optically thick, the thermal 
conditions at the wind affect the hydrostatic structure. This 
is unlike typical stellar systems with optically thin winds. 
As a consequence, the hydrostatic structure has to be solved 
together with the wind, though the governing equations are 
different for the two regions. 



3.2.1 Hydrostatic Region 

In the hydrostatic region, any energy which is generated by 
the viscosity and not advected radially, is transported in 
the vertical direction through either convection or radiative 
transfer, or both. 
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The first equation describing this region is that of hy- 
drostatic equihbrium, 

p dz dz 

The temperature gradient is determined according to 
the energy transfer mechanism. It is given by 

' j-ldPT 



dT 

dz 



= < 



7 dz P' 
' AacT 3 ' 



in the convective zone, 



in the radiative zone, 



(20) 



where p is the density, F is the vertical radiative flux and 7 
is the adiabatic index. Convection is present if the standard 
Schwarzschild criterion is satisfied. But for convection to 
be efficient, the convective flux must be smaller than the 
maximum possible which is given by 
,.3 



- conv,max 



(21) 



and v s is the adiabatic speed of sound. The opacity in the 
radiative zone is taken to be the Thomson opacity, as long as 
the radiative flux is smaller than the critical flux above which 
the atmosphere develops inhomogeneities. As described in 
§2.1, the gas becomes inhomogeneous above the critical flux, 
such that the radiative force exerted on the gas is reduced. 

We assume that the relation between the effective Ed- 
dington factor r e ff = F/FeS and the classical Eddington 
factor T = F/J" E dd is empirically given by 

A 



reft = 1 - 



for r > Tcrit, 



= r for r < r c 



(22) 



r cr ; t is the critical T above which inhomogeneities are ex- 
cited, so the effective opacity for Y > r cr i t is given by 

A 



Keff = «Th 



(23) 



Since we expect a continuous r e fi, A, B and r cr i t satisfy the 
equation r crit = 1— A/Tf lit . From theoretical considerations, 
we take r cr i t ~ 0.8 (Shaviv 2001a). This implies a relation 
between the normalization constant A and the power law B, 
that is given by 



(1 



(24) 



Note that because the behaviour of super-Eddington atmo- 
spheres is expected to depend on only the local conditions, 
the parameters are taken to be location independent. 

3.2.2 Continuum Driven Winds 

As described in §2.2, a continuum driven wind is acceler- 
ated from the region where the density is low enough, such 
that the inhomogeneity based opacity reduction becomes in- 
efficient. In this region, the effective opacity approaches the 
microscopic value, such that the radiative flux is again super- 
Eddington. 

The primary equations describing the wind structure 
are the equation of motion 

dv z _ _dP 
dz dz 

where g z = —d^/dz, and the energy conservation equation 



pv 



P9z 



(25) 



F(z) = F al 



0wind 



GMbh GM B h\ 



R 



(26) 



The index "atm" denote values at the top of the hydrostatic 
atmosphere beneath the wind. Note that we assume here 
that the wind geometry is that of a slab. Namely, we assume 
that z <JC R. This assumption breaks down for high accretion 
rates, at which point the present solution fails. 

From the last two equations, and the assumption that 
k = const., we derive the radiative flux and wind velocity 
as a function of z: 



F{z) = F atm exp 



K^wind (<2atn 



(27) 



GM E 



R, 



1 - 



F(z) 



V R ~ ) 



where 



ind GMbh I Rat 11 



+ f, (28) 



(29) 



m is the ratio between the energy flux needed to accelerate 
the wind out of the gravitational potential well, and the 
radiation flux provided to the wind by the system. 

Another aspect of this thick wind is the location of the 
photosphere. While in the slim and thin disk models the 
photosphere resides where the gas becomes optically thin, 
in our case, the photosphere resides much higher, where the 
wind becomes optically thin. The optical depth of the wind 
is given by 



updz. 



(30) 



Note also that the change in the location of the pho- 
tosphere is accompanied by a decrease in the radiative flux 
emitted from the disk (as energy is used to accelerate the 
wind) , hence, a decrease in the effective temperature. 

As described in §5, the typical ratio we obtain between 
the height of the photosphere and the radius is small, that 
is, z p h(r)/r < 1, as long as the accretion rates are not too 
large, i.e., rh < 20rh cr i t . For higher accretion rates, the wind 
geometry ceases to be disk-like, and the solution described 
here breaks down. 



3.2.3 Photon Tired Winds 

As elaborated upon in §2.3, photon tired winds are formed 
when the available radiative flux at the sonic point is insuf- 
ficient to drive the mass loss driven at the sonic point out of 
the gravitational potential well (see Owocki & Gayley 1997). 
Under such conditions, a layer of shocks forms in which the 
effective sonic point moves upwards, and reduces the actual 
mass loss. Using the results from van Marie et al. (2009), we 
empirically model the maximal </> w ind to be, 



^tiring 



0.2 



( i*Edd ) ' ^ 



(31) 



where (firing =-F/(GM BH /.Ratm). 



3.2.4 Boundary Conditions 



The vertical structure of the disk is determined by the local 
equatorial conditions, i.e., the total pressure P(r, z — 0) 
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Log(r/r g ) 

Figure 2. The mass accretion rate (in units of the critical rate, 
rh/rhcrit) as a function of radius, for different outer accretion 
rates onto a BH with = IOMq. The dotted line denotes 

20m cr it, the solid line denotes I0ih cr it, the short dashed describes 
the 5rh cr it case, and long dashed line accretion with m cr ; t . 



and the density p(r, z = 0). The total radiative flux is a free 
variable. It is determined by the vertical structure which has 
to be solved for given the top boundary condition for the 
radiation field. The latter is the blackbody radiation law, 

F(r,z ) = aT* a . (32) 

Here T c a is the temperature at optical depth r = 2/3, as 
obtained by the outward integration. 



4 NUMERICAL SOLUTION 

The problem we are required to solve is devided into two 
parts, vertical and radial integration. Each radial integration 
step is accompanied by a vertical integration. 

For the vertical integration we take an initial guess for 
the radiative flux and solve eqs. 19-20 for the hydrostatic 
part, and eqs. 25-26 for the wind (equations are integrated 
using the 4 th order Runge-Kutta method), and the validity 
of the boundary condition (eq. 32) is checked. This vertical 
integration is repeated with different values for the radiative 
flux (using the shooting method) until eq. 32 is fulfilled. 

Using the emitted flux F, the mass loss rate 4nr(f> w i n ddr 
and the vertical integrations over the pressure and the den- 
sity, eqs. 7-9, 11 and 13 are solved. This integration is taken 
up to the inner sonic radius (where v r = v s ), which should 
reside between n s and r m b- The location of the sonic point 
is determined by the initial guess for the outer pressure, 
p(r O ut,0). A shooting method is used to determined the 
outer pressure. 




L°g(r/r,) 

Figure 3. The vertical height (z/r) of the disk as a function 
of radius, for the same accretion rates as before. The thick lines 
denote the position of the photosphere when a wind is absent, 
while the thin lines denote the photosphere in the thick wind 
when present, in which case the thick line denotes the location of 
the sonic surface. The line types arc as before. Note that for the 
two lower accretion rate m = rh cr i t and m = hm CT it the wind 
is always optically thin and the photosphere coincides with the 
sonic point of the wind for all radii. 




Log(r/r g ) 

Figure 4. The ratio P=P g /Pt t vs. log(r/r g ). At large radii, 
the solutions approach the standard solution, and the radiation 
pressure is less important. Closer to the black hole, the radiation 
pressure becomes more important and the radiative flux increases. 
Near the inner radial sonic point ds/dr > 0, and the gas pressure 
becomes important again. 
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5 RESULTS 

In this section we present the results for accretion disks with 
a = 0.001 around a stellar BH of mass 10M Q and around 
a super-massive BH of mass 1O 6,5 M . We define the outer 
accretion rate as the rate of mass entering the disk at large 
radii before any mass loss takes place. 

For the stellar BHs, we take outer accretion rates of 1, 
5, 10 and 20 rh cr it. We find that the disks lose considerable 
amounts of mass through a wind, such that the real accretion 
rates (= mass passing the sonic radius per unit time, rh rea i) 
are 0.74, 2.4, 3.8, and 5.7 respectively, i.e, in stellar BHs 
with outer accretion rates rh > 5rh cr i t more than half of the 
mass entering the accretion disk in the outer radius will not 
be accreted into the black hole but leave as a wind. Mass 
accretion rates as a function of the radius are shown in fig. 
2. The vertical heights of the disks are given in fig. 3, where 
thin lines denote the location of the photosphere inside the 
wind, while the thick lines denote the sonic surface (the base 
of the wind). At large radii, where the wind is absent, thick 
lines give the heights of the photosphere as well. Note that 
it is impossible to obtain a situation in which all the mass 
accelerates in a wind since part of the energy is used to heat 
the gas, implying that without any net accretion, there will 
be an insufficient amount of energy to drive the matter back 
to infinity. 

The gas to total pressure ratio, /3 in the equatorial plane 
(z — 0), is shown in fig. 4. Apparently, the gas pressure 
dominates at large radii. However, radiation pressure be- 
comes progressively more important at smaller radii. Radi- 
ation pressure is also dominant for higher accretion rates. 
Close to the radial sonic point the gas pressure becomes 
important again. 



5.0.5 Efficiency of Accretion 

The efficiency of accretion is defined by 



Substituting m rea i into (33) we obtain the accretion efficien- 
cies (see table 1). For low accretion rates, the efficiency is 
slightly higher than the pseudo-Newtonian rjo case because 
the wind is optically thin, while for higher accretion rates 
the wind is optically thick and n/r/o < 1. 

5.0.6 Spectra and Luminosities 

The total luminosities leaving the disks in units of Eddington 
luminosity are 0.8, 2.65, 3.9 and 4.85 respectively. Thus, for 
outer accretion rates of about 1.5 m cr it, the total radiated 
luminosity exceeds the Eddington "limit" . 

The accumulated luminosity L(r) = J r ^ 4mrF(r)dr is 
shown in fig. 5. Assuming a local Planck distribution and 
using the effective temperature T e g(r), the emergent lumi- 
nosity per unit frequency, L„, is 

/Tout 

L„=tt 4nrB„(T cS (r))dr (34) 

where B v is Planck function. This is given in fig. 6. 
Table 1 summarizes the above results. 



5. 0. 7 Sensitivity to Model Parameters 

There are two main uncertainties in our model. The first is in 
the basis of all a-disk models, the a parameter. The only the- 
oretical limit on it is that a < 1. We checked the sensitivity 
of our model using three different values of a, a = 0.001, 0.01 
and 0.1 in an accretion disk having m out = 10rh cr ; t around 
a stellar BH of mass Mbh = IOMq. Increasing the a pa- 
rameter by a factor of 10, decreases the total luminosity and 
the real mass accretion rate by approximately 10%. 

The second uncertainty is the exact structure of the 
porous atmosphere which determines the value of the effec- 
tive opacity. The effective opacity is given by eq. 23, where 
B is a free parameter, and A is determined through eq. 24. 
For the same disk as before, we examined the influence of 
different values of B on the disk structure and luminosity. 
The values of B taken are B = 0.5, B = 1.0 and B = 1.5. In 
this case, increasing B by a factor of 2 (i.e., decreasing n e s), 
causes a small increase of 6% in the real accretion rate and 
a modest increase of 15% in the total luminosity. 

In summary, the theoretical uncertainties in the model 
do not translate into large uncertainties in the model pre- 
dictions. Figs. 7-9 summarize the aforementioned results in 
more detail. 

5.0.8 Super Massive BH 

In addition to accretion disks around a stellar BH, we also 
examined super critical accretion around a super massive 
BH of mass 10 6 ' 5 M Q , with outer accretion rates of Wm cr it 
and 20rh cri t. 

We find that the real accretion rates for these disks 
are slightly smaller when compared to the stellar BH case, 
2.6rh cr it and 3.8rh cr ; t respectively. The relative vertical 
height (z/r) is smaller as well (reaching a maximum value 
of z/r « 0.55 for the m out = 20m crit , as compared to 
z/r « 0.65 in the equivalent stellar BH case. 

The main difference is in the emergent spectra. The 
spectra of disks around massive BHs has a maximum at 
lower frequencies (far-UV compared with the x-ray for stellar 
BHs, see fig. 10 for comparison). The total luminosities are 
2.9LEdd and 3.7LEdd respectively, i.e., modestly less super- 
Eddington. 



6 DISCUSSION 

In the present analysis, we searched numerically for possi- 
ble solutions describing super-critical accretion disks around 
black holes, while allowing for the formation of "porous" 
layers with a reduced opacity. We found solutions with sig- 
nificantly super-critical accretion rates, in which the vertical 
disk height is smaller than the radius, that is, slim disks. Be- 
cause the super-Eddington state excites a strong wind, the 
actual mass accretion onto the BH can be notably smaller. 
Solutions were found to exist with accretion rates ranging 
between about 0.5 rh cr i t to about 20m crit . In all cases, there 
is a photon tired continuum driven wind. 

At the low range, the disks are overall sub-Eddington, 
but locally the flux can surpass the critical value, and there- 
fore it can excite continuum driven winds. Namely, the crit- 
ical accretion rate given by eq. 2 is not the lower limit for 
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riiout /m cr it 


ril r eal/ m crit 


Ltot/^Edd 


Energy Dissipation 


efficiency (77/770) 


1 


0.74 


0.79 


0.8 


1.05 


5 


2.4 


2.65 


2.75 


1.1 


10 


3.8 


3.9 


4.45 


1.03 


20 


5.7 


4.85 


6.9 


0.85 



Table 1. Summary of the slim accretion disks around IOMq with a = 0.001. The mass accretion rates are given in units of the critical 
rate while the luminosity and energy dissipation (which is the total energy produced by viscosity) are given in units of the Eddington 
luminosity of the BH. The last column is the efficiency of the disk compared to the standard disk efficiency (see eq. 2), calculated with 
the real accretion rate. The efficiency calculated using the outer accretion rate is smaller. At high accretion rates, a significant wind is 
accelerated, taking some of the energy and reducing the efficiency. 




2 3 
Log(r/r g ) 

Figure 5. Radially integrated luminosity log(L(r)) vs. Log(r/r g ) 
for the same disks as before. 




16 17 

LogMHz]) 

Figure 6. The total luminosity log(i/L„) vs. log(^), for the same 
disks as before. 




Figure 7. Comparison between the vertical structure, z/r, of 
disks with different atmospheric effective opacities and different 
a, for 10rn cr it around Mbh = IOMq. Solid, dotted and short 
dashed lines denote a = 0.001, 0.01 and 0.1 respectively, while 
B = 0.5. Long dashed and dot short-dashed are for a = 0.001 and 
B = 1.0, 1.5 respectively. Thick lines depict the height of the pho- 
tosphere when an optically thick wind is absent (at outer radii), 
and the sonic point (inner radii). Thin lines mark the height of 
the photosphere when present in a thick wind. 



wind generation. At the low range, almost all the energy dis- 
sipated is either radiated away or transferred to the wind. 
Moreover, the wind is then optically thin and the photo- 
sphere coincide with the sonic point of the wind. 

As the outer mass accretion rate increases, the winds 
become more massive, thereby reducing the fraction of mass 
accreted onto the black hole. Also, some of the dissipated 
energy is then advected with the flow into the BH. 

For very high accretion rates surpassing about 20rh cr it, 
the photosphere which resides in the wind is found to be 
located at z > r. This implies that the 1D+1D type of so- 
lution described here breaks down. Instead, one has to look 
for a solution in which the semi-hydrostatic inflow has a 
disk-like solution, while the super-sonic outgoing wind has 
a spherical-like solution. The description of such extremely 
high accretion rate disks is the subject of a future publica- 
tion. 
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Log(r/r g ) 

Figure 8. Comparison between the accretion rates of disks with 
different atmospheric effective opacity parameters, B, and differ- 
ent viscosity parameters a, as labeled in fig. 7. Note that increas- 
ing a decreases rn rea i while increasing B increases m rea [ , but the 
overall effect is not very large. 




LogMHz]) 

Figure 9. Comparison between the spectra, log(yL„ [erg/ sec]) 
vs. log(v[H z]) , for the different conditions given in fig. 7. 

One of the uncertainties in the model is the opacity law 
behaviour of the porous atmosphere. The reasons it is not 
known well is because it depends on the nonlinear radia- 
tive hydrodynamic configuration the atmosphere will reach, 
and unfortunately, there are still no numerical simulations 
or empirical data which can constrain the effective opacity 
law. ft is for this reason that we parameterized the effective 
opacity (see eq. 23) 

As can be seen in figs. 8-9 almost all disks characteristics 




LogMlIzJ) 

Figure 10. The emitted spectra, \og(vL v [erg /sec}) vs. 
\og(u[Hz\), by accretion disks onto a super-massive BH, with a 
mass of 10 6 5 Mq. The outer accretion rates are 10rh cr it (solid 
line) and 20m cr it (dotted line). The horizontal thin lines are the 
total luminosities and dashed line is Eddington luminosity for 
comparison. 

have either a small or a modest sensitivity to the changes in 
the opacity law. For example m re ai/riiout varies between 0.43 
to 0.44, or the total luminosity varies between ~ 0.570£Edd 
to 0.576£e<m, while changing the opacity parameter B by 
50%. This implies that the uncertainties do not undermine 
the model predictions. On the other hand, it would be im- 
possible to use super-Eddington accretion disks to constrain 
the relevant parameters. 
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